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We propose a compact perturbative approach that reveals the physical origin of the singularity 
occurring in the density dependence of correlation energy: like fermions, elementary bosons have 
a singular correlation energy which comes from the accumulation, through Feynman “bubble” di¬ 
agrams, of the same non-zero momentum transfer excitations from the free particle ground state, 
that is, the Fermi sea for fermions and the Bose-Einstein condensate for bosons. This understanding 
paves the way toward deriving the correlation energy of composite bosons like atomic dimers and 
semiconductor excitons, by suggesting Shiva diagrams that have similarity with Feynman “bubble” 
diagrams, the previous elementary boson approaches, which hide this physics, being difficult to do 
so. 


I. INTRODUCTION 


The energy of N interacting elementary bosons has 
been studied in the late 50’s by Brueckner and Sawadai^, 
and by Lee, Huang and Yang^i^ through a mean-field 
procedure which transforms the two-body Hamiltonian 
into a quadratic operator easy to diagonalize by using a 
Bogoliubov-like transformation. Their most striking re¬ 
sult is that the A^-boson correlation energy is singular 
with a dependence on density n = N/L^ in instead 
of n^. More precisely, the ground-state energy of N in¬ 
teracting bosons reads as 
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where m is the boson mass and a is the scattering length 
of the boson-boson potential at hand. Actually, the 
above result with the scattering length appearing also 
in the term has been obtained by Lee, Huang and 
Yang, but not by Brueckner and Sawada. The reason 
is that the former authors use a pseudo-potential which 
reads in terms of a; so, they do not have to bother 
about generating the scattering length in the correlation 
term. The drawback of these two previous approaches 
is that they rely on a mean-field approximation which 
completely hides the physics of the correlation energy 
singularity. A more elaborate field-theory procedure has 
been later proposed^. It provides a better control of the 
performed approximations, but through a more compli¬ 
cated procedure that still hides the physical origin of the 
singularity. 

The purpose of this paper is to reveal that, in spite 
of the different quantum nature of the particles and the 
different singular dependence in the correlation energy 
of bosons and fermions^, namely and (Inn), the 

physics producing these singular dependences is just the 
same: the accumulation of the same non-zero momentum 
transfer excitations from the free particle ground state, 
that is, the k = 0 boson condensate or the 0 ^ |k| ^ kp 
Fermi sea. This understanding paves the way toward 
deriving the correlation energy of N composite bosons 


made of fermion pairs by suggesting the appropriate set 
of Shiva diagrams^ that have similarity with the Feyn¬ 
man “bubble” diagrams leading to the singular corre¬ 
lation energy of elementary bosons. Composite bosons 
of major present interest are atomic dimers^Ti^ made of 
different species of cold fermionic atomsiSr— , semicon¬ 
ductor excitonsH^— made of electrons and holes, and 
polaritonsi^^— which are linear combination of photons 
and excitons. Cold atoms^i^i and polaritons^^ have been 
used as a testbed for low-energy (Bogoliubov) excita¬ 
tions in the mean-field framework of elementary bosons. 
Whether composite bosons defy this mean-field descrip¬ 
tion because of their composite nature remains under de¬ 
bate in spite of the fact that, because of the fermion in- 
distinguishability, one cannot construct an effective po¬ 
tential between composite bosons that is valid beyond 
first order in interaction^. As a direct consequence, el¬ 
ementary boson approaches based on the existence of a 
boson-boson potential cannot be duplicated for compos¬ 
ite boson systems. 

To show the analogy between the correlation energies 
of elementary bosons and elementary fermions, we pro¬ 
pose a compact perturbative approach to the energy of 
N quantum particles that allows catching the physics 
of the singular interaction processes. Once selected and 
summed up, these singular processes lead to the /V-boson 
energy given in Eq. ©• Unlike previous methods, the 
perturbative approach that we here propose can be di¬ 
rectly extended to composite bosons which interact not 
only through the fermion-fermion interactions between 
their fermionic components, but also through fermion ex¬ 
changes. 

The paper is organized as follow: 

• We first give some general arguments for understand¬ 
ing what should be and what really is the density depen¬ 
dence of the correlation energies for N elementary bosons 
and N elementary fermions, in order to understand why 
these density dependences end up being different, albeit 
produced by the same physical processes. 

• Next, we propose a compact perturbative approach 
to derive the A^-boson energy, which allows catching the 
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physics of its various terms in a transparent way. We also 
provide the key commutators which enable us to calculate 
these terms easily. 

• We then recover the Wboson energy given in Eq. (ED 
through the explicit summations of the ladder diagrams 
associated with the scattering length, and of the bubble 
diagrams associated with the correlation energy singular¬ 
ity. We also discuss the required cancellation of overex- 
tensive contributions that come from disconnected dia¬ 
grams, as standard in perturbative expansion. 

• We conclude with the state-of-the-art for compos¬ 
ite boson systems and the fundamental problem raised 
by the Pauli exclusion principle between the fermionic 
components of composite quantum particles. 

In the supplemental material, we outline the original 
derivations of the A'^-boson energy proposed by Brueckner 
and Sawada, and by Lee, Huang and Yang, the physics 
responsible for the correlation energy singularity being 
hard to catch from these mean-field approaches. The 
supplemental material also contains some heavy diagram¬ 
matic parts associated with high-order perturbative ex¬ 
pansions, which support our procedure but are not nec¬ 
essary to follow its spirit. 

II. GENERAL ARGUMENTS 

We look for the solution of 

{H-£nMn) = 0, ( 2 ) 

where H = Hq -|- U is the Hamiltonian of the quantum 
particles at hand. The free part is given by 

Ho = '^ EkC^Ck (3) 

k 


the extra 2 coming from spin. By contrast, since (6^)^ ^ 
0, two bosons can be in the same state. So, the iLg ground 
state corresponds to the so-called “condensate” with all 
bosons in the lowest energy state. 

To understand that the difference in the density de¬ 
pendence oi £n for bosons and fermions comes from the 
same physics but within different Hg ground states, let us 
expand Sat in terms of the number of particles involved 
in scattering processes, namely 

£n = {--- )N+{- • • )iV(iV-I)+(- • • )N{N-l){N-2)+- ■ ■ . 

( 7 ) 

The fact that the Fermi sea extension depends on fermion 
number through Eq. © brings density-dependent pref¬ 
actors into the above expansion, while such case does not 
exist for bosons. So, the resulting density dependences of 
£n/N for bosons and fermions have to be different, even 
if they are induced by same physical processes. 

This is already seen from the N term which corre¬ 
sponds to the energy in the absence of interaction. The 
energy of N free bosons in the k = 0 condensate reduces 
to 0 while the energy of N fermions in the Fermi sea 
scales as Nep oc 

To analyze the other terms, we use dimensional argu¬ 
ments. As the Wparticle energy is an extensive quantity, 
£n /N only depends on density. Potential scatterings de¬ 
pend on sample volume as while each sum over 

momentum brings a factor when transformed into an 
integral. Let us successively consider what these com¬ 
ments imply on the boson and fermion energies. 

A. Bosons 

• The N{N — 1) term of £is[ comes from (1, 2, 3, • • •) 
interactions between two bosons taken among N. 

As these two bosons come from the condensate, process 
involving a single potential scattering can only appear for 
q = 0. So, it gives, within a numerical prefactor. 


with Ek = k^/ 2 m while the interaction part can be writ¬ 
ten as 

^ j 4+q4'-qCk'Ck . (4) 

q kk^ 

The elementary particle operators fulfill 
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with 77 = —I for bosons, Ck = bl, and 77 = -|-1 for 
fermions, cl = al- Since (al)^ = 0, two fermions cannot 
be in the same state; so, the Hg ground state corresponds 
to a Fermi sea having all k states with 0 ^ |k| ^ kp oc¬ 
cupied, the number of 3D fermions with up and down 
spins being related to the Fermi momentum kp through 


7V = 2 ^ 1 = 

O^k^kF 


37r^ 


( 6 ) 


N{N — l)vo oc Nn. ( 8 ) 

Processes involving two potential scatterings must con¬ 
tain a q sum to have the proper extensivity, and an en¬ 
ergy denominator to be an energy-like quantity. This 
energy denominator can only be the energy of the boson 
pair (q, — q) excited from the condensate, as a result of 
the q scattering. This gives, within a numerical prefac¬ 
tor, 

v'^ 

7V(7V-l)y]y^ oc Afn. (9) 

Processes involving three potential scatterings must 
contain two q sums, and two energy denominators. And 
so on. 

This shows that the Nn terms of £n come from inter¬ 
action between two bosons excited out of the condensate 
through the so-called “ladder processes”, each additional 
Vq being accompanied with a q sum and a 2eq denomina¬ 
tor, in order to make this additional part dimensionless 
and sample volume free. 
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Such ladder processes physically correspond to the 
renormalization of the Vq scattering appearing in Eq. 
through the repeated excitations of one boson pair from 
the condensate, this pair changing from (qi,—qi) to 
(q 2 , —q 2 ), and so on. These ladder processes change the 
bare first-order term of Eq. m into 

iV(iV - l)wo cx iVn , (10) 


the renormalized q == 0 scattering being commonly writ¬ 
ten in terms of the scattering length a as 


47ra 


( 11 ) 


• The 7V(7V — l)(iV —2) term of £n comes from interac¬ 
tion between three bosons. To be an energy-like quantity 
with proper extensivity, this term must contain to ^ 2 
interactions, (to — 1 ) energy denominators, and (to — 2 ) 
sums over momentum. 

The term with two scatterings does not exist, because, 
in the absence of q sum, it can only contain the q = 0 
process; so, the required energy denominator 2eq would 
be zero. 

The term with three potential scatterings, which must 
have one q sum and two energy denominators, reads, 
within a numerical prefactor, as 


iV(iV-l)(fV-2)^-^cxfVn 


( 12 ) 


The A^-boson correlation energy results from interac¬ 
tion between more than two bosons. It fundamentally 
reads, within a renormalization of the potential scatter¬ 
ing through ladder processes, as 


qi 


(16) 
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We then note that, for r’q-j.o ^ 0, the above qi sums 
diverge in the small qi limit, with similar divergences 
occurring for more than four bosons. As standard, the 
summation of these singular terms overcomes their diver¬ 
gences and produces a finite contribution to the A^-boson 
energy with a density dependence somewhat larger than 
its Nn'^ first term. To show that the summation of these 
singular processes leads to a dependence requires 
the knowledge of the numerical prefactors in Eq. (HU. 
These prefactors, which depend on the number of ways 
these scattering processes appear, are obtained from a 
precise counting that cannot be done from dimensional 
arguments only. 

We conclude from the above analysis that the repeated 
interaction between two bosons taken from the A-boson 
condensate produces a dressing of their bare interac¬ 
tion through a set of ladder processes like the one^^ of 
Fig.m the resulting linear term in density being given by 
Eq. (HOl). 


The term with four potential scatterings, which must 
have two q sums and three energy denominators, reads 


A(A- 1)(A- 2) 
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(13) 

Although not obvious from dimensional arguments alone, 
these four potential scatterings split as written above, the 
q 2 sum belonging to the ladder processes that renormal¬ 
ize the Uq interaction between two bosons according to 
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(14) 


Higher-order terms in this potential renormalization 
come from more than four scatterings between two 
bosons. 

• The A(A— 1)(A — 2)(A— 3) term of £n comes from 
interaction between four bosons. The energy-like contri¬ 
bution with proper extensivity involves four interactions 
at least. It reads 

iV(A-I)(A-2)(A-3)E7:^«lV^'- (15) 

qi 

More interactions between four bosons lead to a renor¬ 
malization of the potential scattering, as in the case of 
three bosons. And so on. 



FIG. 1: (a) Ladder process between two bosons from the con¬ 
densate. (b) Feynman diagram obtained by closing the 0 
lines. 

Higher-order terms in density come from interactions 
between more than two bosons. The most singular ones, 
given in Eq. (HSl), correspond to repeatedly exciting and 
de-exciting the same (qi, — qi) pair from the condensate 
through “bubble” diagrams like the ones shown in Fig. [21 

B. Fermions 

• The A(A — I) term of £n comes from interaction 
between two fermions taken from the Fermi sea. 

The energy-like term involving a single interaction 
scales as 

A(A-I)ukp OC , (17) 

for Vq taken as the 3D Coulomb scattering, Uq = 
for q yf 0, and vq = 0, the q = 0 process 
being eliminated for electrons in a positive ion jellium 
which ensures the system neutrality^. Since only q ^ 0 
processes exist, this single-interaction term comes from 
fermion exchange inside the Fermi sea. 
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FIG. 2: Bubble processes resulting from three interactions 
between three bosons 0 excited from the condensate. They 
provide equal contributions to the A'^-boson correlation en¬ 
ergy. 
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FIG. 3: (a) Direct Coulomb process and (b) exchange 

Goulomb process between two fermions from the Fermi sea. 


The energy-like term involving two interactions must 
have one energy denominator which scales as i ^^tl 
one q sum which brings a N factor. So, it scales as 

7V(iV-l)|iV^| oc Nn°. (18) 

With two interactions, there in fact exist direct and ex¬ 
change processes, as shown in Fig. O The direct process, 
associated with a two-bubble Feynman diagram, is sin¬ 
gular due to the infiniteness of its two Uq scatterings in 
the small q limit. 

Ladder processes similar to the ones existing between 
two bosons like Fig. [T] also exists between two fermions. 
However, as each interaction comes with a Uq scattering 
which scales as Ukp: energy denominator which scales 

as Ekp : and a q sum which leads to a factor N, each 
rung of a ladder process scales as 

oc (19) 

which is small in the large density limit. This is why lad¬ 
der processes and scattering length are never considered 
for Coulomb interaction between electrons in the dense 
limit, while they play a key role for bosons in the dilute 
limit. Note that this is not so for electron-hole systems 
due to the poles these ladder processes produce that are 
associated with exciton formation in the dilute regime. 


• The N{N — 1)(A^ — 2) term oi £n comes from in¬ 
teraction between three fermions. The term with two 
potential scatterings does not exist for the same reason 
as that for bosons. The next-order term with three po¬ 
tential scatterings, two energy denominators, and one q 
sum, scales as 

iV(iV-l)(fV-2)|iV^| oc Nn-^/^. (20) 

The energy-like term involving four fermions, four scat¬ 
terings, three energy denominators, and one q sum, scales 
as 

fV(iV-l)(iV-2)(iV-3)|iV^l oc 7Vn-2/3,(21) 

I ^kF J 

and so on. 

Again, direct and exchange processes exist for terms 
involving more than two fermions. Precise calculations^ 
show that, much like in the case of bosons, the numer¬ 
ical factors of the direct terms diverge in the small mo¬ 
mentum transfer limit. When summed up, these diver¬ 
gences transform the nP dependence of the first term, 
given in Eq. dl), into a somewhat larger (In n) depen¬ 
dence. These singular terms correspond to bubble pro¬ 
cesses, just like the singular terms for bosons that leads 
to a dependence for the correlation energy. This 
understanding is commonly known for fermions but not 
so for bosons. The unique but major difference between 
fermions and bosons is that, for fermions, these singular 
terms are associated with excitations from a Fermi sea, 
whereas, for bosons, the singular terms are associated 
with excitations from a condensate with all particles in 
the k = 0 ground state. 

Before turning to the procedure we propose to derive 
the Wboson energy given in Eq. (ED. we wish to make an 
important comment. In the above arguments, we start 
with the free quantum particle ground state and we per¬ 
form similar expansion in terms of the number of particles 
involved in scattering processes. Yet, we end up with a 
small density expansion in the case of bosons, and with 
a large density expansion in the case of fermions. The 
reason is that, in order for the free particle ground state 
to be a good starting point, this state must not be very 
much changed by interaction. Such is the case for a very 
large Fermi sea, i.e., a very dense fermion system. By 
contrast, the Wboson condensate is extremely narrow 
in energy by construction; so, to have it not very much 
changed by interaction, the number of scattering events, 
I.e., the boson density, must be very small. 


III. EQUATION FULFILLED BY THE A-BOSON 
GROUND-STATE ENERGY 

We want to solve Eq. © in the case of bosons. For 
N/L^{= n) and V small enough, the H ground state 
stays close to the Hq ground state, that is the A-boson 
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condensate IOat) = where \v) denotes the vac¬ 

uum state. We force this IOa?) state into the problem by 
inserting 


|0Af)(0jv| 

(OatIOat) 


( 22 ) 


in front of I'i/’Af) bi Eq. ([2]) and we multiply the resulting 
equation either by IOat) or by P±. As Ho\0n) = 0, we 
get 


0 = (OatIE - £n\0n) + {0n\VP^\M ■ (23) 

\UAr|t)Ar/ 

As P±Ho\On) = 0, multiplication by P± yields 


where s refers to the perturbative order in the W inter¬ 
action. So, = £^\ except for s = 1. 

Equation (|M1) readily gives 0 = = £^'^ and 

0=A^^; so, 

^(1) N{N-1) (O^lEolOw) 

Because W operators only enter Eq. (1^ . higher-order 
terms in interaction come from q 0 processes. 

As A AT is second order at least in W, the two next- 
order terms read as 

A^^\0n\0n) = {0n\WPj_^P^W\0n) , (31) 

—Ido 


0 = P^HoP±\M+P±V\0N)^^f^ 

+P±{V - £n)P±\'iPn) , (24) 


from which we get 


P±\M 


P±^P±V\0n) 

— Hq 


(OAflV^Af) 

(OatIOat) 


+P±^P±{V - £n)P±\M , 

—Ho 


(25) 


that we iterate. By inserting the resulting P±\'tpN) into 
Eq. (1^ and by noting that Pj_£at|0at) = 0, this gives 
the equation fulfilled by £1 at as 


00 / 1 AS 

0 = {0N\{V-£N)Y.iP±^P±{V-£N)] IOat). 

(26) 

To solve the above equation in an easy way, we split V 
as Vb + 14^ where 

^=EYE^U^kT-q^k'6k. (27) 

kk' 


A^^(OatIOat) = {0n\WP^^P^WP^^P^W\0n) , 

— ilQ —-tlQ 

(32) 

while, for s > 4, they have a more complicated form; 
Aj^^(0AT|0Y = {0n\WP^^P^/^^Pa_^P^W\0n) , 

—lio —-tlQ 

(33) 

the first operators being given by 

4"^ = WP^^P^W - (34) 

—Ho 

jP = WP^^Pi^WP^-^PiW (35) 

— Uq —tlQ 

-aP {p^P^^p^w + WP^^^P^ - aP . 

From the above two equations, we see that disconnected 
processes exist when more than three IE’s act inside 
the IOat) condensate, as standard for perturbative expan¬ 
sion. They generate overextensive contributions, which 
are eventually canceled out by the A at parts of as 
shown more in details below. 


contains all non-zero scattering processes. As Vo act¬ 
ing on any state having N bosons gives voN{N — l)/2, 
Eq. (I^Sl) ultimately appears as 


00 

0 = (0^|(1E-A^)^ 

s=0 



1 


P±{W-An) 


where An = £n — voN{N — l)/2. 


IOat) , 
(28) 


V. SOME USEFUL COMMUTATORS 


(s) 

The A]y quantities are easy to calculate with the help 
of the following commutator 



q#o 


(36) 


IV. POTENTIAL EXPANSION 


where Tq = Yk^k+q^k- This excitation operator is such 
that 


We solve Eq. ((28l) as a IE expansion with A at written 
as 


T 51 
^ q? 


= b 


t 

p-^q ’ 


An = 




s=0 


(29) 


from which we readily get 

TqlOjv) = A'fc^lO^-i). 


(37) 


(38) 
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By iteration, the above equations lead to 


VF|0^) 




IOat-i) 


iV(iV- 1) 
2 


q#o 


(39) 


where creates a finite-momentum pair from 

the condensate. Interaction involving such a pair follows 
from 


As previously shown (see Eq. (l30l) '). the A^-boson en¬ 
ergy at first order in interaction is equal to VoN{N — l)/2, 
while all higher-order terms involve ^ 0 processes 
through W interactions. 

• The second-order term, given by 


£ 


( 2 ) 

N 


^(2) N{N-1) 

— ^5 


E 

qi#o 


~2£qj^ 


oc Nn , 


(42) 


is visualized by the diagram of Fig. [S] 



q#0 

(40) 


VI. STRUCTURE OF THE POTENTIAL 
EXPANSION 



The calculation of Ajv requires the knowledge of s op¬ 
erators W acting on the [Oat) condensate, (OatIVE'^IOat) 
containing the parts of IF'*[Oat) in which all excited 
bosons are back into the condensate: 

The first operator W excites two q = 0 bosons from 
the condensate to become a (qi,—qi) pair, as shown in 
Fig. m Equation (l3^ readily gives 


FIG. 5: Two W interactions inside the IOa^) condensate re¬ 
duce to the ladder process of diagram (a). Diagram (b) is 
the standard bubble diagram for the same process, this Feyn¬ 
man diagram being obtained by closing the “open” lines 0 of 
diagram (a). This process appears with a N{N — l)/2 prefac¬ 
tor due to the number of ways to choose the two interacting 
bosons 0 among N. 


P±^P±W\0n) 

—Ho 


N{N-1) 

2 



(41) 

The N{N— 1)/2 prefactor comes from the ways to choose 
the two interacting bosons 0 among N. 





0 


N(N-\) 


0 


FIG. 4: One W acting on IOa^) brings one boson pair out of 
the condensate since W contains qi yf 0 processes only. There 
are N{N — l)/2 ways to choose the two interacting bosons 0 
among N. 

A second W operator acting on this excited state can 
do three things: 

(i) It can induce an additional q 2 scattering between 
the same two bosons, leading to terms in i?qj^_|_q 2 | 0 Ar- 2 ) 
with the same N{N — l)/2 prefactor. 

(ii) It can induce an interaction with a third bo¬ 
son 0 from the condensate, leading to terms in 
fc^lqi-qjfcijifcijjOAr-a) with a prefactor [V(V-1)/2](V- 
2), from the {N — 2) ways to choose the third boson 0. 

(iii) It can excite a second boson pair (q 2 , —q 2 ) from 
the condensate, leading to terms in B^^B^JOn- 4 ) with 
a prefactor [N{N - 1)/2][(7V - 2){N - 3)/2]. 

More W interactions bring more bosons 0 from the 
condensate into play. Some examples of the resulting 
FE^IOat) states are given in the supplemental material. 


( 3 ) 

• The third-order term, , has two origins: 

(i) One contribution comes from the ladder process 
shown in Fig. |6]in which (qi, q 2 , qa) are different from 0 
but qi -b q 2 + qa = 0. For qa set equal to q 2 — qi, this 
contribution reads 


N{N-1) 

2 


E E 

qiEo q^Efo.qi) 


( 2eq')( 2£qj) 


OC Nn. (43) 


0 

0 





0 

0 


(a) 


FIG. 6: Ladder process resulting from three W interactions 
between two bosons 0 from the IOac) condensate. 



(ii) The second contribution comes from the bubble 
processes involving three bosons 0, shown in Fig. [2] It 
reads 


N(N-l), 

2^--^(iV-2) 






oc Nn‘‘ 


(44) 


qiEo 


the extra factor 2 coming from the two processes in this 
figure which give equal contributions. 

When compared to the second-order term shown in 
Fig. El we see that a third interaction FE either adds 






























7 


one more interaction between two bubbles as in Fig. |6]or 
adds a third bubble, as in Fig. [51 The former term, which 
appears with a N{N — l)/2 prefactor, corresponds to a 
renormalization of the interaction between two bosons 
through ladder processes. The latter term, which appears 
with a prefactor N{N — 1){N — 2), is part of the correla¬ 
tion energy. This term actually diverges in the small qi 
limit for 0 because it contains the 3D integral 

of (l/qf)^. The Nn? dependence of this singular three- 
bubble process is transformed into a dependence, 

when joined with similar bubble terms of higher order in 
W. To understand how this happens, we must consider 
higher-order terms in W. These higher-order terms more¬ 
over bring a new feature: the appearance of overextensive 
contributions coming from disconnected processes. The 
various fourth-order terms in W are given in the supple¬ 
mental material, as well as the cancellation in the large 
JV limit, of their overextensive parts. To explicitly show 
the necessary cancellation of these overextensive parts 
through the Ajy part of becomes more and more 
cumbersome when s increases. From now on, we focus 
on extensive contributions resulting from connected dia¬ 
grams. 


VII. LADDER PROCESSES 


At any order in W, there are terms that come from 
the repeated interaction between two bosons, as shown 
in Figs. EllSl and [H The sum of these ladder processes 
leads to the effective scattering hq shown in Fig. [71 which 
is the solution of the integral equation 


Va = Va 


^a — c 


qVO 




-Va' 


(45) 


This summation is commonly associated with the two- 
body scattering length a defined as 

lim hq = 47r . (46) 

q-^o ^ mL^ 


It is possible to rewrite the energy term coming from 
any number of interaction between two bosons taken from 
the condensate, that is, all processes like the ones of Figs. 
|5l ini and [T] plus the first-order term given in Eq. (1501) , in 
terms of the scattering length a, as 


iV(iV-l) 27r 

---'Co = N—na , 

2 m 


(47) 


in agreement with the first term of Eq. ©• 

A similar interaction renormalization exists between 
bubbles. It leads to replacing Cq with Vq. Actually, this 
simple procedure is only valid in the small q limit, which 
is when bubble processes are most singular; for larger 
q’s, the interaction renormalization is affected by the sur¬ 
rounding excited states. 

To understand this important point, let us consider 
three-bubble processes. The one involving three Cq’s, 



FIG. 7: Effective scattering resulting from the repeated inter¬ 
action between two bosons. 


shown in Eig. [51 leads to the contribution given by 
Eq. dill). These three Cq interactions are also renormal¬ 
ized by adding more W interactions, as shown in Fig. jS] 
The three bubble diagram then transforms up to fourth 
order in W into 



^iV(iV-l)(iV-2)^^^|\ (48) 


The Vq interaction associated with Fig.jSKc) does not read 
exactly the same as the Vq interactions coming from the 
diagrams in Fig. |5Ka,b). However, the dominant contri¬ 
bution of these three-bubble processes comes from small 
q’s because of the divergence of the l/(—2eq)^ factor it 
contains; for such momenta, Vq ~ Vq. 


0 0 



FIG. 8: Repeated interaction between three bubbles leading 
to Eq. (l48ll . 































VIII. ORIGIN OF THE SINGULARITY 


We now consider the V-boson correlation energy, i.e., 
the terms coming from interactions between more than 
two bosons, and we explain why the dominant correlation 
energy term does not depend on density as n^, but as 
. This non-perturbative result originates from small- 
q divergences such as those in Eqs. (jH) and (l44ll . The 
standard procedure to remove these singularities is to 
sum up similar divergent terms. 

In order to explicitly show how this summation re¬ 
moves the singularities, let us consider the lowest-order 
terms represented by the diagrams of Figs. [51 [H and 21 
of the supplemental material. In this set of processes, 
the system has one pair (q, — q) excited from the k = 0 
condensate, along all intermediate steps, this pair being 
ultimately de-excited back into the condensate. These 
diagrams lead to 


— 


iV(iV- 1) 




q/O 




q#o 


(-2eq)^ 


+4iffi^(]V-2)(iV-3)y: 


q#0 






(Nv^y 




(49a) 

(49b) 


for N large. The effect of summing up this series is clear: 
it brings an interaction term to the pair kinetic energy, 
which stays finite in the small q limit, thus serving as a 
small-q energy cut-off, provided that vq is finite. 

Of course, the above result follows from appropriate 
prefactors. Let us show that bubble processes involving 
s ^ 3 interactions lead to 






(50) 

Since this term involves s bosons, its N factor is easy 
to understand: N{N — l)/2 comes from the number of 
ways to choose the hrst pair of bosons 0 from N identical 
bosons 0 in the condensate, as shown in Fig. 01 The 
other N factors come from adding additional bosons 0 to 
this initial pair, the number of bosons that remain in the 
condensate decreasing as {N — 2), (iV — 3), ... 

The 2®“^ prefactor is more tricky. It comes from topo¬ 
logically different diagrams, such as the two diagrams in 
Fig. m that involve s = 3 interactions. Let us provide a 
systematic way to generate this 2^~‘^ prefactor. We start 
with the q interaction between two boson lines, shown in 
Fig. m and we introduce two operators U and V : the U 
operator adds an up-line connected to the upper boson 
line by a q interaction, as in Fig. Hla); the V operator 
adds a down-line connected to the lower boson line by a 
q interaction, as in Fig.jHlb). To construct the diagrams 
having s bubbles, we apply {U + 'Dy~‘^ to the two-boson 




(a) 


9 


<1 


(b) 


FIG. 9: The U operator adds an up-line to the diagram of 
Fig- SI as in (a), while the T) operator adds a down-line, as in 
(b). 


diagram of Fig. [Hand we close the uppermost and lowest 
boson lines by a —q interaction. This binomial expansion 
leads to a prefactor 2®“^, since each term in the expan¬ 
sion corresponds to topologically different process, while 
the integral is the same. The corresponding third-order 
and fourth-order bubble diagrams are respectively shown 
in Fig. 121 and 21. 

By rescaling Cq in Fq. (I49b|) as 2Nvqx, we get 


A)v — 


1 

2 



4TrmV2rn{2Nvoy^^SN , 


(51) 


where, for Vq taken equal to vo when 0 < Eq < and 
0 otherwise, in order to insure the large-q convergence 
of the q sums which is normally insured by the natural 
decrease of Vq, the Sn sum appears as 


Sjv 


p^/‘ 2 Nvq 

/ dx 

lo 

\fx 


2a: -|- 1 


1 n 

I 

1 n 

V 2Nvo 

^f2 \ 

-1 

o 


(52) 


The first term of 5'^?, which produces a Vq term in A^y, 
comes from the first term in Eq. (I49al) . It should be 
subtracted off because it is already included in the ladder 
series that produce the part of Ajv linear in density. By 
letting the potential cut-off ft go to infinity in the second 
term of S'at, we get a contribution in to A^y, that 
can be written in terms of the scattering length at first 
order in interaction, oi ~ mL^vo/A'K, as 

ivl^(na?)3/2. (53) 

maf 

Difference of this result from the second term of the N- 
boson energy in Eq. (HD is two-fold: (i) it reads in terms 
of the first-order scattering length m, instead of the full 
scattering length a; (ii) its numerical factor is different. 

(i) To get the full scattering length a, instead of oi, in the 
above equation is quite easy: we just have to note that 
ladder processes also exist between bubbles, as shown in 
Fig. [S] In the small q limit where bubble contributions 
are singular, they lead to replacing vo in Eqs. (ED and 
(l52)l by Vo, that is replacing oi by a in Eq. ()53p . 

(ii) The numerical factor in Eq. ()53p differs from the one 
of the term in Eq. (HD, by 4% onlyi. This evidences 
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that the dominant contribution to the A^-boson correla¬ 
tion energy comes from the repeated excitation of just 
one boson pair from the condensate. The missing small 
contribution comes from processes in which more than 
one boson pair (q, — q) are excited from the condensate. 
Actually, the correlation energy appears as 


A 


corr 

N 


1 1 I 1 i 

2 ^0 -2^'! - 2^^q 2 (-2£q - 2Nv^r 



(54a) 

2^^^^(—2f —27V?; Us-i 

(54b) 


• The first term of Eq. (I54al) comes from the excitation 
of one boson pair, as given by Eq. (|49bp . 

• The second term comes from the process shown in 
Fig.iini in which two boson pairs are excited and then 
de-excited, the energy denominator —2£q being replaced 
by —2eq— 2Nvq, to account for the renormalization from 
bubble processes, as shown in Eq. (I49al) . Since at least 
one excited boson pair must remain in the intermediate 
steps, there is only one way to destroy the two pairs: they 
are created first and then destroyed. 



FIG. 10: Process involving the excitation of two boson pairs 

(q,-q)- 

• The third term of Eq. (l54aD comes from the two pro¬ 
cesses shown in Fig. 1111 in which three boson pairs are 
excited and then de-excited. These two processes corre¬ 
spond to the two ways to destroy three boson pairs: (i) 
one can excite three bosons pairs and then destroy them 
all, as in Fig. [TlT al: (ii) one can excite two boson pairs, 
destroy a pair, excite another pair, and finally destroy 
them all, as in Fig, [m b'). 

Following this simple rule, it is easy to count all possi¬ 
ble ways of destroying s excited pairs. The fs prefactor 
in Eq. (I54b|) represents the number of ways that s excited 
boson pairs can be destroyed. Brueckner and Sawadai 
have shown that this fs factor is just the numerical factor 
of the generating function 

OO 

(55) 

S = 1 

Using this function, it is possible to show that the energy 
given in Eq. (1541) is nothing but the correlation energy 


0^ 

0^ 


0^ 

-q;, 

0^ 


0^ 

0^ 




(a) 


-^0 

-^0 


-^0 

^0 


r^O 

s? 

^0 


0^ 

0^ 


0^ 

-q; 

0^ 


0^ 

0^ 



^0 

^0 

^0 

^0 

^0 

^0 


(b) 


FIG. 11: Processes involving the excitation of three boson 
pairs (q, —q). Process (a) excites three pairs and then destroy 
them all. Process (b) first excites two pairs, then destroys one 
pair, then excites another pair, and finally destroys them all. 


obtained from the mean-field approach, in this way re¬ 
covering the precise numerical prefactor for the term 
given in Eq. (HD- 


IX. STATE-OF-THE-ART FOR COMPOSITE 
BOSONS 

The major problem when dealing with composite 
bosons (“cobosons”) comes from the fact that, due to 
possible fermion exchanges, there is no way to assign 
a given fermion pair to a composite boson. As a di¬ 
rect consequence, there is no way to write an effective 
Hamiltonian with an interaction potential between co¬ 
bosons that is valid beyond first-order processes, in spite 
of what all bosonization procedures intend to do. So, 
previous approaches used for elementary bosons, which 
are based on the existence of a boson-boson potential 
such as the one of Eq. ®, cannot be used for composite 
boson systems. Of course, it is always possible to stay 
with th e po tential between elementary fermions, as done 
in refs. However, the ladder processes that lead to 
the formation of a composite boson, have to be selected 
and summed up first, before considering interaction be¬ 
tween composite objects. In this regard, the coboson 
many-body formalisni^ proposed a decade ago is far sim¬ 
pler because it avoids the handling of ladder processes 
associated with composite boson formation. 

The dominant way^ two cobosons interact is through 
the Pauli scattering induced by the Pauli exclusion prin¬ 
ciple, as represented by the Shiva diagram in Fig. &). 
Since, in it, no interaction takes place, this dimension¬ 
less scattering is by construction missed by all bosonized 
Hamiltonians, because in them, only enter energy-like 
quantities. 

Attempts have been made to reach the iV-coboson en¬ 
ergy linear in density in the case of cold atom dimers 
within a contact potential^i^, and in the case of semi¬ 
conductor excitons2ii2^ within the long-range Coulomb 
potential. The cold atom problem is rather simple be¬ 
cause in the linear term in density, only enter two bosonic 
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'>c" 

(a) 


r—777^1:1 

-<- f 

(d) 

FIG. 12: (a) The Pauli scattering A(^^). (b) The direct 

Coulomb scattering (c) The “in” exchange Coulomb 

scattering (d) The “out” exchange Coulomb scatter¬ 

ing r\Z)- 



i n -. 


-t^.i 


j m 




(b) 


atoms; so, this linear term can be obtained by numer¬ 
ically solving a four-fermion problem. The next-order 
density term, that is, the correlation energy, is far more 
tricky to get because it requires solving a full many-body 
problem. Indeed, it is reasonable to expect that, as in the 
case of elementary bosons, this term has a singular den¬ 
sity dependence that comes from more than three fermion 
pairs. So, a numerical brute-force calculation is hopeless. 

In the case of cold atomic gases with a contact potential, 
it has been showit^, using a procedure similar to the one 
proposed by Lee-Huang-Yang, that the A^-dimer energy 
is also given by Eq. o, where a is the scattering length 
between two dimers. 

These procedures, hard to extend to the long-range 
Coulomb potential, fail to address the main physical 
problem when dealing with composite bosons: how do 
fermion exchanges enter the result? The repeated in¬ 
teraction between two cobosons must contain the direct 
Coulomb scattering shown in Fig. [T^ bb as well as the 
two exchange Coulomb scatterings shown in Fig. lT^ c.db 
It can also contain the Pauli scattering induced by 
fermion exchange, multiplied by the difference of the co¬ 
boson energies of the process at hand in order to get an 
energy-like quantity. So, we expect the effective scatter¬ 
ing appearing in the ladder processes between two co¬ 
bosons to read2^ 

c{^;) = + + 

+ (• ■■){Em + Er,-E,- E,)\ (™ *) (56) 

It is far from obvious that the set of exchange processes 
between two cobosons will appear just the same in the 
case of many cobosons. One simple reason is that Pauli 

( m i \ 

ra i J or more also 

exist^. So, even if the next-order term in density still be¬ 
haves as its prefactor may not solely depend on 

the scattering length between just two cobosons. The 
understanding, through the present work, that the 
singular dependence comes from the accumulation of the 
same non-zero momentum transfer excitations from the 


condensate provides a great clue and a guidance for sug¬ 
gesting the series of Shiva diagrams that should lead to 
the most singular term in the A^-composite boson cor¬ 
relation energy, and for possibly confirming the result 
obtained in the case of atomic dimers^. 


X. CONCLUSION 

We propose a compact perturbative procedure to de¬ 
rive the energy of N interacting bosons obtained by 
Brueckner and Sawade^i and by Lee, Huang, and Yang^i^, 
through the summation of an infinite series of diagrams 
of the “ladder” type to get the scattering length and 
of the “bubble” type to get the correlation energy — 
which depends singularly on the density n as while 
its density dependence is (In n) for fermions. One ad¬ 
vantage of the present procedure is that it provides a 
unified way to derive the singular correlation energy of 
fermions and bosons in terms of q ^ 0 excitations from 
their corresponding non-interacting ground states: the 
Bose-Einstein condensate for bosons and the Fermi sea 
for fermions. The other advantage is that it gives a phys¬ 
ical picture for how the singular behavior arises from the 
accumulation of q ^ 0 excitations. This understanding 
will serve as a valuable guide to study composite boson 
systems such as Bose-Einstein condensates of cold atoms 
or semiconductor excitons, by selecting appropriate Shiva 
diagrams having similarity with “bubble” Feynman dia¬ 
grams. 
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XI. SUPPLEMENTAL MATERIAL 


Appendix A. Brueckner-Sawada approach 

Brueckner and Sawada proposed a mean-field approach 
to the energy of N interacting elementary bosons. Their 
boson Hamiltonian reads E[ = EIq -f V with Hq given in 
Eq. (3) and V given in Eq. (4), the operator cj^ being the 
elementary boson creation operator They performed 
a mean-field treatment in which the creation and destruc¬ 
tion operators 5 q and bo are replaced by a scalar v^Aq, 
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with A^o being the number of bosons 0 in the condensate. 
The total number of bosons then reads 

iV = iVo + ^iVk = iVo + ^6],6k. (A.l) 

k5.^0 kT^O 

In the N ~ Nq regime, the V dominant terms have 
four and two boson 0 operators, but not three due to 
momentum conservation. The part in vq reads as 

^ k^O k^O ^ 

(A.2) 

due to Eq. (|A.1I) . while the part in Uq^o is given by 



Since ek=o = 0, we end with a mean-field Hamiltonian 
given by 


The simplest way to estimate this energy is to force 
the scattering length into the problem. It reads, up to 
second order in interaction. 


47r 


a 

niL^ 


Vo ^ no 


+ E 

q#o 


V- 


‘’-2eq '*■ 


(A.IO) 


This allows us to write the ground state energy as 

= N^na + E'j^ . (A.ll) 

with n = N/ and E'j^ given by 

= 2 E (\/^q - (A-12) 

n^n “I 

^ 4:TTmV2m{2Nvo)^^^K. 

For Vq taken equal to vo up to q infinite and for Sq 
rescaled as 2NvqX, the K integral reads as 

AT = y dx\/x(^x'^ + X — {x + -) + —^ . (A.13) 


L 

27r 


= ^N^ + ^J2hq, (A.4) 

q#0 


where hq appears as 

hq = eq(^bl^bq -I- &Lq6_q^ -|- t'ql^&q&Lq + bqb-qj . (A.5) 


By noting that the primitive of this integral is given by 

/fM = l{(*+l)V“-xV“}-|{(x+l)’/“+ix»/^}+ A , 

(A.14) 

so that K{0) = —4/15 and K{oo) = 0, we ultimately find 
the A^-boson energy in the mean-held approximation as 


The effect of the k = 0 condensate is to renormalize 
the q yf 0 interaction as z/q = Nvq and to dress the 
q 7 ^ 0 boson energy as Eq = Eq -I- zzq, the interaction 
ultimately appearing as creation or destruction of (q, —q) 
boson pairs. 

The hq Hamiltonian is easy to diagonalize in terms of 
Bogoliubov-like operators 

/3q = ^q cosh 9q + b-q Sinh 6q . (A.6) 

For tanh20q = z/q/sq, we hnd that hq can be written as 


p(MP’) 

N 


27r 
—n 
m 



(A.15) 


This result differs from Eq. (1) through the fact that the 
term depends on the scattering length oi associated 
with the bare potential vo, instead of a associated with 
the interaction scattering dressed by ladder processes. 
Approaches beyond the mean-held approximation are re¬ 
quired to transform ai into the full scattering length a in 
the correlation term. 


~ cosh 20q ) + COShEq + ^-q^-q) ' 

(A.7) 

So, the mean-held Hamiltonian ultimately reads 

= 4“^) + ^ ^e~ 2 - Z.2 /3t, (A.8) 

q#0 


where /8^ creates excitation with energy y Eq — t'q from 

the mean-held ground state. The mean-held ground state 
energy is given by 




q#o 


^q “ ^q “ ^q 


(A.9) 


Appendix B. Lee-Huang-Yang approach 

The Lee-Huang-Yang approach^i^ takes the force be¬ 
tween two bosons as a hard-sphere interaction. This 
hard-sphere interaction is further approximated by a 
pseudo-potential21 


y(pseudo) ^ 4 ^ ^ 


i>j 


dr. 


(B.l) 




where is the position of the i boson and r^- is the dis¬ 
tance between i and j bosons. For such a hard-sphere 
interaction between two bosons with small momenta k 
(|k|a <C 1), the parameter a physically corresponds to the 
hard sphere diameter. When the potential in Eq. dBH) 
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is used in the case of just two bosons, one finds that a 
corresponds to the full two-body scattering length. So, 
the great advantage of the above pseudo-potential is to 
force the full scattering length into the problem from the 
very beginning. Besides the delta function (5(ri — r^) de¬ 
scribing a contact force, the peculiar form of {d/drij)rij 
reflects the boundary condition of hard spheres at r = a, 
below which the wave function must reduce to zero; it 
guarantees that there is no l/r singularity when r —>■ 0. 
When written in second quantization, the potential in 
Eq. (IB.II) appears as 


state gives, using Eqs. (39) and (40), 

WBl\0^_,) = ^ (C.2) 

q2 7^o 

+ (iV-2) ^ 

q2 7^o 

q2#o 

This state is shown in Fig. [T31 It contains states having 


= Sf Si T/ s i: . 


k'k q 


The Hamiltonian of hard-sphere bosons is then treated 
in a mean-field approximation like in the Brueckner- 
Sawada approach. By singling out the dominant inter¬ 
actions involving four 0 bosons and two 0 bosons, the 
approximated Hamiltonian reduces to 


H 


- -p^ -I- - Inn —r N 

mL^ 2 r-^o Or 

qT^o 


X 



47ra7V\ 

mL^ J 


{b\bq + 6)lq6_q) 


AiraN 

mL^ 





(B.3) 


?l 

9i+ 92^-|-f^ 

^2^ -q^ 

-9|-92^^- < ^ 

(a) 


0 

0 


N(N-l) 

2 


92 

9i 

-9,-92 



92^ 
—92 —^ 
9,^ 
-9.^ 


0 92 —-0 (^^-2) 

924 9i 

0 9, - 92 < ^ — <ry — 0 N(N - 1) 

0 - 9 ,^-®— 0 

(b) 


92 



0 (jV-2)(jV-3) 

0 2 

0 N(N-l) 

0 2 


FIG. 13: A second W acting on the state of Fig. 4 can act on 
the excited pair, as in (a), or can involve one or two additional 
bosons 0 from the condensate, as in (b) and (c). 


which can be exactly diagonalized by using a standard 
Bogoliubov procedure. We then find that the ground- 
state energy reads as Eq. (|A.12p . with fq replaced by 
Eq + AiraN/mL^, and t'q replaced by AiTaN/mL^. In 
addition to having the full scattering length a appearing 
into the problem at no cost, another advantage with the 
pseudo-potential (IB.Ill is that the last term of Eq. (IA.12|1 , 
which now reads (47raiV/mL^)2/2eq, appears in a natural 
way to remove the 1 /q^ singularity that is hidden in the 
first term. 

Appendix C. IE® acting on IOat) 

• A first IE potential acting on [0^,) excites one boson 
pair from the condensate (see Eq. (39)) 

w\0^) = I0A2-2), (C.l) 

qi#o 

where = bJj^bLq^ creates a boson pair. The resulting 
state is shown in Fig. 4. 

• A second IE acting on this one-boson-pair excited 


two, three and four bosons excited from the condensate, 
with prefactors N{N — I)/2, [N{N — I)/2](A — 2), and 
[A(A—1)/2][(A —2)(A—3)/2], respectively, which come 
from the number of ways to choose these bosons 0 among 

N. 

• We now consider a third IE acting on the states of 
Fig. [ 13 ] having two, three, and four excited bosons, and 
we successively consider what IE does on these states. 

(i) IE acting on one excited pair 

Adding a qa 7 ^ 0 interaction to the diagram of 
Fig- USKa) is just the same as adding a q 2 7 ^: 0 interac¬ 
tion to the diagram of Fig. 4: we simply have to replace 
the qi scattering by a series of (qi,q 2 ) scatterings. As 
(qi,q 2 ,q 3 ) differ from 0, the process in Fig. [Hl al is the 

(3) 

only ladder process that contributes to A)^'. 

(ii) IE acting on three excited bosons 

Adding a qa 7 ^: 0 interaction to one of the two diagrams 
in Fig.lWbl. which are topologically equivalent, leads to 
the diagrams in Figs. (US and (HU). 

— In Fig. [151 the qa interaction takes place between 
two of the three bosons already involved in Fig. lT^ bb So, 
these processes also appear with a [N{N — 1)/2](A — 2) 
prefactor. As (qi,q 2 ,q 3 ) differ from 0, the processes of 
Fig.[IKa,c) do not contribute to A^^\ By contrast, since 
we can take qi = q 2 = —qa, the process of Fig. [TSI bl 
contributes to 
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FIG. 14: A third qa 7 ^ 0 added to the processes in Fig. [13] 
can involve the same two bosons 0 as in (a), a third boson 0 
as in (b), or excite another boson pair from the condensate 
as in (c). 
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FIG. 15: The qs interaction takes place between two of the 
three bosons involved in Fig. Illl bl. 


FIG. 16: The qa interaction involves one or two additional 
bosons from the condensate. 
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— In the processes of Fig. [I 6 i:a,b,c), the qa 7 ^ 0 inter¬ 
action involves a fourth boson taken from the condensate, 
while in Fig. [TBr dl. this interaction excites a (qa,—qa) 
boson pair. So, the former processes appear with a pref¬ 
actor N{N — 1){N — 2){N — 3), while the latter process 
appears with a prefactor N{N — 1){N — 2){N — 3){N — 4). 
However, since qa ^ 0, all these processes do not con¬ 
tribute to . 

(iii) W acting on two excited pairs 

Finally, we can also add one qa 7 ^ 0 interaction to the 
two excited pairs shown in Fig. HHc). 

— This can be done inside a pair, as in Fig. HZIa), 
or between two pairs, as in Fig. I17f bl. However, since 
q 2 7 ^ 0 , these processes do not contribute to A}^ . 

— The third W interaction can also involve one more 
boson from the condensate, as in Fig. I17l cl. or involve 
two more bosons, as in Fig. fTTT dl. However, these pro- 
cesses do not contribute to A]^ in the same way, because 
(qi, q 2 , qa) differ from 0 . 


FIG. 17: Third-order interactions involving four or more 
bosons. 


Appendix D. Fourth-order term in W 

The fourth-order term in W is obtained by adding 
a q 4 7 ^ 0 interaction to the third-order state FF^IOat), 
shown in Appendix C. Processes giving a non-zero contri¬ 
bution to (Oat IFF^ I Oat) involve two, three or four bosons 
0, while those involving five bosons reduce to zero, be¬ 
cause FF only contains momentum transfers that differ 
from zero. 


1. Two bosons involved 

We can add a q 4 7 ^ 0 interaction to the ladder process 
in Fig. 6 . For q 4 = —qi — q 2 — qa, this process gives a 
non-zero contribution in the IOat) condensate, as shown 
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in Fig. 1. By setting q2 = q .2 ~ ‘li Qs = Qs ~ ^ 2 ’ its 
contribution to reads as 

N{N- 1) ^ 

qi5^0 q^5^(0,qi) 

(xNn. (D-1) 

This term contains four fq’s, which bring four 1/L^ fac¬ 
tors, three q sums, which bring three factors, and a 
N{N—l)/2 prefactor; so it leads to a contribution in An. 
Together with the term shown in Fig. 6, this part of 
participates in the renormalization of the Vq interaction. 


E 


^-q3^q3-q2^q2-qi^' 


qi 


qs^C^iq^) 


(-2eqO(-2£qA(-2£qJ 



FIG. 18: This process, when combined with the three-bubble 
diagram in Fig. 2(a), participates in the renormalization of the 
Wq interaction between bubbles through the ladder processes. 


0 




(b) 


FIG. 19: More complicated four-interaction processes be¬ 
tween three bubbles. 


contains four t^q’s, two q sums, and a N{N — 1){N — 2) 
prefactor; so, it leads to a contribution in Nn^. 

We can also have more complicated three-bubble pro¬ 
cesses, in which the added interaction links different bub¬ 
bles, as shown in Fig. ITW a). Finally, we can add two 
interactions and one bubble to the two-bubble diagram 
of Fig. 5, as shown in Fig. lTW b). The processes shown in 
Fig. fTOl a) and fTOl b) actually contribute to higher-order 
density terms^i^. 


3. Four bosons involved 

With four W interactions involving four bosons 0, we 
can have the four-bubble diagram as shown in Fig. 



FIG. 20: Four-bubble diagrams: the system stays with one 
excited boson pair (qi, —qi) across all the intermediate steps. 

These bubble processes produce a contribution to 





(b) 



(d) 


£. Three bosons involved 


FIG. 21: Topologically different diagrams corresponding to 
the four-bubble diagrams in Fig. 1201 


We can add one interaction to the three-bubble pro¬ 
cesses shown in Fig. 2 by repeating one of the interac¬ 
tions, as in Fig. 1181 So, it participates in the renormal¬ 
ization of the Uq interaction between bubbles. This term 


which reads as 




^-qAqi 2 

oc Nn . 


qi#o 


(-2SqJ^ 
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the extra 4 factor coming from the four topologically dif¬ 
ferent processes shown in Fig. \TI\ 



FIG. 22: Disconnected processes resulting from four W inter¬ 
actions between four bosons 0 from the condensate. 

We can also have four bosons involved through the 
disconnected processes shown in Fig. They lead to 


These two terms contain four Uq’s, two q sums, and a 
N prefactor coming from the number of ways to choose 
the two boson pairs (0, 0) among N. So, they lead to 
an overextensive contribution in , which is canceled 
out by the part of in Eq. (34). Let us show it. 

Appendix E. Cancellation of overextensive terms 

Up to now, we have shown that the 7V-boson energy 
at first order in interaction is equal to vqN{N — l)/2 
(see Eq. (30)), while all higher-order terms depend on 
q ^ 0 processes through the W interaction. given 

in Eq. (42), has a unique contribution that scales as Nn, 
while has a contribution, given in Eq. (43), which 
also scales as JVn, and another contribution, given in 
Eq. (44), which scales as Nn^. In the same way, the part 
of coming from the first term of in Eq. (34) 
has extensive contributions in JVn (Eq. (|D.1I) 1 and in 
Nn^ (Eq. (ID.21) 1. However, the part of also has two 
overextensive contributions, given in Eq. (ID.31) . Actually, 
these two terms can be combined into 


N(N- 1) (A - 2)(A- 3) 


2 2 

q^^oq'^^o (“2eq) (-2eq - 2eq/)(-2eq) 


+ EE 75 


V-aV. 


qt^q' 


q/OqVO 


( 2£q') ( 2eq 2eq/) ( 2£:q) 


N(N- 1) (A-2)(A-3) 


'^q ^q' 

^ f-2e„12 2.- (-2eq0 ' 


qT^O ^ q'7^0 


(D.3) 

oc . 


(E.l) 

It then becomes easy to see that, when inserted into 
Eq. (34), this overextensive part is canceled out in the 
large N limit by the part of 


'N ■ 
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